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Abstract 

We consider a non-interacting one-dimensional gas accelerated by a constant and 
uniform external field. The energy absorbed from the field is transferred via elastic 
collisions to a bath of scattering obstacles. At gas-obstacle encounters the particles 
of the gas acquire a fixed kinetic energy. The approach to the resulting stationary 
state is studied within the Boltzmann kinetic theory. It is shown that the long 
time behavior is governed by the hydrodynamic mode of diffusion superposed on a 
convective flow. The diffusion coefficient is analytically calculated for any value of 
the field showing a minimum at intermediate field intensities. It is checked that the 
properly generalized Green-Kubo formula applies in the non-equilibrium stationary 
state. 
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1 Introduction 



The problem of constructing adequate statistical ensembles for non-equilibrium 
stationary states (hereafter denoted by NESS) has a long history (see e.g. 
[1,2]). Related to it there are very interesting findings concerning the nature 
of fluctuations in the presence of spatial gradients inducing stationary flows 
[3]. The computation of exact expressions for the probability distribution of 
simple NESS allows to study different static properties as, for example, the 
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long range density correlations [4]. The recent development of the theory of 
fluidized granular matter has also attracted a lot of attention to the structure 
of NESS (see for example [5,6,7]). Finally, there is at present an important 
research carried on in order to properly define the entropy production accom- 
panying stationary dissipative currents (see e.g. [8]). 



However, whereas the modes of approach of fluids to thermal equilibrium 
are well understood, there is still no satisfactory theory of the dynamics of 
reaching NESS. In both cases the final state is stationary, but in NESS, in 
contradistinction to equilibrium, there persist non-vanishing fluxes sustained 
by the coupling to environment. 



The kinetic theory of fluids revealed the fundamental role of the separation 
between different time scales in the approach to equilibrium (see e.g. the pre- 
sentation of Bogoliubov's ideas in [9]). In particular, in the evolution of the 
density of particles in the one-particle phase space all the degrees of freedom 
but those represented by the densities of globally conserved quantities relax 
at a short time scale producing a state close to local equilibrium. The subse- 
quent slow evolution is essentially governed by changes of the hydrodynamic 
fields (an original construction of the hydrodynamic modes within the kinetic 
theory has been developed by P. Resibois [10]). One can wonder whether this 
kind of mechanism based on the time-scale separation still persists when the 
final state involves some dissipative stationary process related to the current 
flowing through the system. This is precisely the question we want to address 
in this paper. 



Our object here is to exploit an analytically soluble one-dimensional model 
discussed in [11] in order to perform a detailed study of the evolution towards a 
NESS. Although the system is very simple indeed, we find that the possibility 
of a rigorous study of its dynamics is precious from the point of view of the 
search for the laws governing the approach to NESS. The coupling to the 
environment is represented in the model by the action of an external field. One 
of the original results of the paper is an explicit formula for the field dependent 
coefficient of diffusion governing the evolution of the only hydrodynamic mode 
in the approach to the stationary state. 



In Sec. 2 we define the system and its dynamics. The following analysis of the 
approach to the stationary state is presented in Sec. 3. Finally, Sec. 4 contains 
the concluding comments. 
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2 One-dimensional system and its dynamics 



We consider non-interacting point particles of mass m moving in one dimension 
with acceleration a under the influence of an external constant and uniform 
field. The particles are surrounded by a bath of point obstacles of mass m 
that move with velocities ±£7, both directions being equally probable. At 
elastic collisions between the particles and the obstacles a simple interchange 
of velocities is taking place. The scattering obstacles are uniformly distributed 
with a number density n. 

The statistical state of the particles is described by the probability density 
f(r, v, t) for finding a particle at time t at point r with velocity v. We assume 
that the evolution of the density f(r,v; t) is governed by the linear Boltzmann 
equation [11] 

(d d d \ f 

dt + V dr + a dv J ^ V) ^ = n J d ^ V " °' ^ C ' ~ ^ V ' *^ c ^ ' 

(1) 

putting 

4>{v) = \ [5(v + U) + 5(v - U)] (2) 

for the velocity distribution of the obstacles. When writing equation (1) it has 
been supposed that the particles always encountered obstacles with velocities 
distributed according to (2). The possibility of recollisions with obstacles per- 
turbed by previous collisions with the particles is not taken into account by 
the Boltzmann collision term in (1). 

In dimensionless variables w = v/U, x = nr, and r = Unt equation (1) takes 
the form 

^[5{w - l)/i-(x, t) + S(w + 1)(JL+(X, t)], (3) 

where F(x,w,t) = f(x/n,Uw,r/nU)U/n is the dimensionless probability 
density, and 

A t ±(^, r )= J dw\w ± l\F(x,w,r). (4) 

In equation (3) there appears the dimensionless intensity of the field e = 
a/nU 2 . 

Provided that appropriate boundary conditions are supplied at ±oo, equa- 
tion (3) has a non-equilibrium stationary solution representing a homogeneous 
NESS whose velocity distribution F (w; e), that can be chosen to be normal- 
ized as / dwF (w;e) = 1, has been derived and analyzed in [11]. The NESS 
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is characterized by a constant current Vness(c) = / dwwF (w; e) that shows a 
linear response, Vness( 6 ) ~ e, in the weak field limit |e| <C 1 and a non-linear 
response, Vness( c ) ~ sgn(e) -1/ 1 e | , for strong fields |e| I. 



Whereas the paper [11] focused on the properties of the stationary distribu- 
tion, our aim here is to analyze the approach to the asymptotic homogeneous 
NESS starting from an arbitrary inhomogeneous initial condition. This re- 
quires solving the time-dependent Boltzmann equation. The way towards the 
construction of the solution has been already found in [1 1] . Using the Fourier- 
Laplace transformation 



F(k,w,z) = J™ dr J dxe- ZT ~ lkx F{x,w,T) 
one finds from (3) the integral equation 



(5) 



F(k, w, z) = H(k, w, z) + — — J- exp 



X(w) _ ik w2 ~ 1 



x | exp 
+ exp 



z(w 



l + z(w + l) 



2e 

8[e(w - l)]jl-(k,z) 
6[e(w + l)]/t+(M) 



where 



x (w) = [( w + l)\w + l\ + (w- l)\w - 1|] /4, 
that obey the relation x{ w ) — ~x{~ w )i 



(6) 
(7) 



fi±{k,z) 



dw\w ± l\F(k, w, z), 



and 



_ poo 

H(k,w,z)= e- ZT - ikx H(x,w,r) (9) 
J o 

relates to the initial condition F(x, w, 0) through 

H(x, w, t) = exp[(x(w — er) — x{w))/e]F(x — wr + er 2 /2, w — er, 0) (10) 

Eq. (6) is implicit expressing F(k,w,z) in terms of fl±(k,z) that are linear 
functionals of F(k,w,z). However, the solution can be made explicit in a 
straightforward way. Indeed, multiplying Eq. (6) by \w ± 1| and integrating 
with respect to w one obtains a closed system of linear equations for jl± in 
the form 



jl-(k,z) 




h_{k,z) 






h+(k,z)_ 



Here 



M(k,z; e) 

z (k, z) — J dw\w ± l\H(k,w, z) 



h 



(11) 
(12) 
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and M(/c, z; e) is a two by two matrix analytic in z given by 



M(k,z;e) 



M 11 (k,z;e), M 12 (k,z;e) 
M 21 (k,z; e), M 22 (k,z; e) ; 



(13) 



with 

Mn = 

M 12 : 

M 21 = 
M 22 -- 

When e > 0, the matrix M(fc, z; e) takes the form 
M(A;,z;e) 



— — 7 1 dw\w — 1 


exp 


xH 

e 


, w 2 - ll 

— «K 

2e 




exp 


xH 

e 


., w 2 -l 

— ik 

2e 


— : — - / dwroi 1 


exp 


xH 

e 


., w 2 -l 

— ik 

2e 


— — 7 ( dw\w + 1 


exp 


xH 

e 


. lW 2 -l 

— ik 

2e 



exp 
exp 
exp 
exp 



1 — z(w - 


1)1 


e 




1 + z(w 


+ 1) 


e 




1 — z(w - 


1)1 


e 




1 + z(w 


+ 1) 


e 



0[e(«; - 1)] 

0[e(iu + l)] 
0[e(«;-l)] 
0[e(iu + l)]. 



1 - {A - 5)e( 2z+1 )/ 2e , C-D + (B- A)e-( 22+3 )/ 2e 
-(A + 5)e( 22+1 )/ 2e , 1 - {A + £)e-( 22+3 )/ 2e — c — D 



where 

A(k, z; e) 



— / dww exp { 

2e Ji I 2e 



; (1 +ik) + 2wz-ik j 



I f 00 f 1 r 

z;e) = — dw exp < — — w 2 (l + ik) + 2wz — ik 
2e Ji I 2e L 



— / dw w exp { 

2e 7-i F I 2e 



C(/c, z; e) 
D(M;e) = i£^exp{-l 



i£;w 2 + 2(tw + l)(z + l)-ik 
ikw 2 + 2(w + l)(z+l) - ik 



(14) 

(15) 
(16) 
(17) 
(18) 



The corresponding formulae for e < follow from the important symmetry 
relations 



M u (k,z;e) = M 22 (-k,z;-e) 
M 12 {k,z;e) = M 21 (-k,z;-e) 



(19) 
(20) 



3 Modes of Approach to the Non-Equilibrium Steady State 



The linear equation (11) contains all information about the modes by which 
the system approaches the NESS. Whether these modes behave for long wave- 
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lengths as hydrodynamic modes is the question we want to study now. To this 
end we use the explicit solution of the Boltzmann equation. The structure of 
(6) and (11) implies that F(k,w,z) is an analytic function in the complex 
plane except at points z where the inverse matrix M _1 does not exist. Thus 
the zeros of the determinant Det(M) define the singularities in the z-plane 
of the Fourier-Laplace transform F(k,w,z). In particular, isolated zeros Zi(k) 
corresponding to simple poles will produce after applying the inverse Laplace 
transformation an exponential time dependence of the form 

F(k,w,T) = Y / ^e z * (k)T , (21) 



with coefficients on for each mode depending on the initial condition (F de- 
notes the spatial Fourier transform of F). 

We performed a systematic numerical survey of the zeros of the determi- 
nant Det(M) for different values of k and e. It has been found that re- 
gardless of the value of e there is a single isolated zero which behaves like 
z = —Dk 2 — ikV + 0(k 3 ) for k <C 1 . All the other zeroes, both for finite 
k and in the limit k — > 0, have negative real parts and are located outside a 
band of a certain finite width around the imaginary axis. That is, these other 
zeros do not have an accumulation point with vanishing real part. Therefore, 
for long wavelengths, there is exactly one slow diffusive mode (coefficient D) 
combined with convective transport with velocity V . All the others are fast 
kinetic modes. The presence of an unique slow mode is related to the fact that 
the mass is the only conserved quantity. No zeroes with positive real part were 
found, reflecting the fact that the NESS is stable. 

The slow mode can be obtained analytically by inserting the asymptotic for- 
mula z = —D(e)k 2 — ikV(e) into the equation Det(M) = and then solving 
for D and V keeping only dominant terms in k. The relations (20) imply that 
the zeros of the determinant of the matrix M(fc, z; e) coincide with the zeros of 
the determinant of the matrix M(— k, z; — e). It follows the symmetry relations 



D(e)=D(-e)=V(\e\), (22) 
V(e) = -V(-e)= S gn(e)V(\e\). (23) 

One finds then a unique solution 



V(A 4l + e-^] + [e-l-(l + e)e-^]I(e) 
U e[l-e- 2 A] + [-e + 3 + (l + e )e- 2 A]/( e ) 1 > 
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and 



V{e) = e 2/e ( - 2e 3 + e 4 + 2e 5 + 8e 2 /(e) + 2e 3 /(e) - 10e 4 /(e) - 8e 5 /(e) - 10e/(e) 

- lle 2 /(e) 2 - 2e 3 /(e) 2 + 7e 4 /(e) 2 + 6e 5 /(e) 2 + 4/(e) 3 + 8e/(e) 3 + 12e 2 /(e) 3 
+ 14e 3 /(e) 3 + 6e 4 /(e) 3 ) + e 4/e ( - 4e 3 + 2e 4 - 4e 5 + 20e 2 /(e) - 8e 3 /(e) 

- 12e 4 /(e) + 16e 5 /(e) - 20e/(e) 2 - 2e 2 /(e) 2 + 20e 3 /(e) 2 + 22e 4 /(e) 2 - 12e 5 / 
+ 4/(e) 3 + 16e 2 /(e) 3 + 8e 3 /(e) 3 - 12e 4 /(e) 3 ) + e 6 ^( - 2e 3 - 3e 4 + 2e 5 



+ 12e 2 /(e) - 26e 3 /(e) + 22e 4 /(e) - 8e 5 /(e) - 10e/(e) 2 - 3e 2 /(e) 2 + 30e 3 /(e) 2 
- 29e 4 /(e) 2 + 6e 5 /(e) 2 - 8e/(e) 3 + 20e 2 /(e) 3 - 22e 3 /(e) 3 + 6e 4 /(e) 3 ) 



e(e - /(e) - e/(e) + e 2/e (-e - 37(e) + e/(e)) )' 



(25) 



with 



(26) 



The drift velocity V(e) (24) coincides, as expected, with the stationary average 
velocity Vness found in [11]. In Fig. 1 we plotted its dependence on e, showing 
the transition from a linear to a nonlinear response. 



The explicit expression for the diffusion coefficient (25) is quite involved. A 
plot of it is presented in Fig. 1. It should be noticed that D(e) has an interesting 
structure as a function of the intensity of the external field showing a minimum 
for e ~ 1.60. The asymptotic formulae for small and large intensities are 



. ^ 1 15 2 



23(e) 



4 - 7T 

TIP' 



9tt-20 _ 

+ (2tt) 3 / 2 6 



|e| < 1, 
e > 1. 



(27) 
(28) 



It should be remarked that in the limit of vanishing external field the value 
-D(O) coincides with that following at equilibrium (e = 0) from the Green- 
Kubo formula. For large values of e, the velocity of the obstacles becomes 
negligible compared to that acquired by accelerated motion and their action 
resembles that of stopping centers. Therefore by purely dimensional analysis 
one can predict the exponent 1/2 in the dependence of D(e) in (28). 

The appearance of a minimum in -D(e) shown in Fig. 1 is an interesting point 
which calls for interpretation. The minimum can be understood by analyzing 
the partial collision frequencies with obstacles moving with velocities +1 or 
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Fig. 1. Top: Average velocity Vness(c) in the NESS as a function of the imposed 
acceleration e. Bottom: Field dependent diffusion coefficient D (e). 

— 1. In the NESS, these collision frequencies are given by 

A*o±(e) = / dw\w ± l\F (w; e), (29) 

where the average is computed with the use of the NESS distribution F (w; e) 
given in Eq.(40) of Ref. [11]. As shown in Fig. 2, the collision frequency /t - 
with obstacles moving to the right (in the direction of the accelerating field) 
has a minimum at e = 1.37. In order to understand this fact consider a particle 
that has just collided with an obstacle with velocity —1 getting instantaneously 
its velocity. If the field is weak, it will continue moving to the left encountering 
obstacles with velocity +1 and will be sent again towards those with velocity 

— 1 the collision frequency /io- remaining high. Also if the field is large enough, 
the particle will turn rapidly to the right and gain a large velocity leading again 
to a large value of /io-- There are however intermediate values of acceleration 
where the particle remains a long time turning to the right under the action 
of the field before the next collision occurs without getting large velocities and 
thus reducing the collision rate. This fact explains the origin of the minimum 
in /io-(e). Finally, the presence of the minimum in /io- can be related to the 
observed minimum of D(e) because the lowering of the collision rate for a 
range of values of e produces an evolution that is closer to a deterministic one, 
and thus less diffusive. 

Knowing the analytic solution of the initial value problem for the kinetic 
equation (3) we could also verify if the diffusion coefficient (25) follows from the 
Green-Kubo formula when the NESS distribution is used in the evaluation of 



8 



2.4- 



2.2 



2.0 



1.2- 



1.0 



1.4 



1.8- 




0.8, 



'() 



1 



2 



3 



4 



5 



Fig. 2. Partial collision frequencies in the NESS, solid line no + , dashed line no-. 
the autocorrelation function. We thus considered the Green-Kubo expression 



where F GK (w,r) is the solution of the equation (3) with initial condition 
F GK (w, 0) = (w-V)F (w) and V is the NESS drift velocity. The computation 
process is as follows. Given the initial condition F GK (w,0), the equation (11) 
is solved for ji± (note that z must be nonzero because otherwise the inverse 
matrix M(A; = O,^)" 1 does not exist). Then, Fqk{w,z) is obtained and the 
coefficient -Dgk can be found from the relation 



We performed the calculation obtaining a finite value which coincides with 
(25). We thus conclude that when the system is in the NESS, both in the linear 
and in the nonlinear response regimes, the diffusion coefficient is given by the 
Green-Kubo formula provided the NESS distribution is used in the averaging 
process and the fluctuations of the velocity around the non-zero mean value 
Vness are considered. This conclusion is also confirmed by an analogous result 
obtained for a simpler one-dimensional case where only nonlinear transport is 
present [12]. 

Finally, we have checked that beyond the linear response regime the system 
does not obey the fluctuation-dissipation theorem, relating the diffusion coef- 
ficient D and the mobility The mobility is defined as [i = dV^ESs/dt- Eval- 
uating the non-equilibrium stationary temperature T = ((w — Vness) 2 ) using 
the stationary distribution F , it is directly checked that D(e) ^ T(e)/i(e). The 




(30) 



(31) 




(32) 



9 



equality is only satisfied at the equilibrium case e = 0. A similar phenomenon 
has been predicted in granular gases [13]. 



4 Concluding comments 

We have considered a one-dimensional system of particles absorbing energy 
from a constant and uniform external field, and dissipating it via collisions 
with a bath of scattering obstacles with a fixed velocity distribution (2). In 
this situation the appropriate linear Boltzmann equation (3) predicts the ap- 
proach to a non-equilibrium steady state (NESS) characterized by a mass 
flux that can present both linear and nonlinear response depending on the 
strength of the field. It has been found that the evolution of arbitrary initial 
conditions toward the NESS involves short time scale kinetic modes produc- 
ing the stationary velocity distribution followed by the slow long wavelength 
hydrodynamic diffusion superposed on a convective current flowing with the 
average velocity of the NESS. 

Analytic calculations show an interesting structure in the dependence of the 
diffusion coefficient on the intensity of the external field, with a minimum 
at intermediate values of the field. The minimum is related to the presence 
of another minimum in the collision frequency between the particles and the 
obstacles, thus reducing the number of randomizing scattering processes. Fur- 
thermore, we have also shown that the value of the diffusion coefficient could 
be obtained from a Green-Kubo formula by considering the time displaced 
peculiar velocity correlation function (subtracting the average NESS velocity) 
averaged over the NESS distribution. 

It thus appears that the modes of approach to the NESS when only the mass 
remains as a conserved quantity involve the corresponding single hydrody- 
namic mode which is the classical process of spatial diffusion. The diffusion 
coefficient does depend on the external field intensity and can be obtained at 
any value of the field from an appropriately generalized Green-Kubo formula. 
We conclude that the approach to NESS shows qualitative analogy to that of 
reaching equilibrium. Hopefully this conclusion remains valid for a large class 
of systems in which the stationary state results from the balance between the 
energy flow absorbed from an external field and collisional coupling to a kind 
of thermostat. 

The authors thank P. Cordero for helpful discussions. This work has been 
partly financed by Fondecyt research grants 1030993 and 7040123 and Fondap 
grant 11980002. 
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